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Abstract. Regime of MOND, which is used in astronomy to describe the gravitating 
systems of island type without the need to postulate the existence of a hypothetical 
dark matter, is generalized to the case of homogeneous distribution of usual matter 
by introducing a linear dependence of the critical acceleration on the size of region 
under consideration. We show that such the extrapolation of MOND in cosmology is 
consistent with both the observed dependence of brightness on the redshift for type la 
supernovae and the parameters of large-scale structure of Universe in the evolution, 
that is determined by the presence of a cosmological constant, the ordinary matter of 
baryons and electrons as well as the photon and neutrino radiation without any dark 
matter. 
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1. Introduction 

In the framework of modified Newtonian dynamics (MOND) pLj we get quite the 
reasonable theoretically and successful phenomenologically explanation for the empirical 
Tully-Fisher law, that regularly relates the visible masses of spiral galaxies to 
asymptotically constant rotation-velocities of stars in the region of dominating "dark 
matter halo" (see figure [1]) as well as for profiles of rotational velocities of stars 
in galaxies, wherein the regime of flattening the rotation curves does not take 
place. Moreover, the MOND is able to explain many other phenomena caused by 
the inhomogeneity of baryonic matter-distribution in the Universe at the scales of 
gravitationally coupled systems (see reviews in [2l[3lll]). For instance, in disk galaxies, 
wherein a dominant contribution to the baryonic mass is given by the interstellar gas, 
the matter density is determined with a high accuracy, since an uncertainty, caused by 
the extraction of stellar masses from their visible magnitudes, is significantly reduced. 
The MOND gives a precise agreement of its predictions with the experimental data on 
the rotation speeds of matter around the centers of galaxies without any free parameters, 
while the uncertainties of results are caused by errors of measurements, only [5]. 

By construction, the basic conclusion of MOND is the statement that the 
observational data can be confidently described without any introduction of non- 
baryonic dark matter at the galactic scales of matter-inhomogeneity as well as at the 
scales of gravitationally coupled galactic clusters [3] due to the specially fitted universal 
modification of the gravity law if the acceleration of free falling g is less than the critical 




Figure 1. The correlation of visible baryonic masses in the disk galaxies with the 
rotation velocities of stars in the region of flat rotation curves, as shown in comparison 
with the MOND prediction: v'^ = GMg^. The mass is given in units of solar mass 
Mq. The figure is taken from [2]. 
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value Qq ~ 1.2 ■ 10~^°m/s^, so that = GMg^/r'^, where M is the matter mass, while 
r is a distance to a reference point. 

However, such the conclusion conflicts with cosmological modeling the Universe 
evolution and properties of its large scale structure, that is characterized by a high 
degree of homogeneity for the spatial distribution of matter {6p/p ~ 6T/T ~ 10~^ 
for fluctuations of energy density p and temperature of cosmic microwave background 
radiation (CMBR), T), because such calculations lead to the necessary introduction of 
dark matter with the density being approximately 5 times greater than the baryonic 
density [6], if we follow the equations of general relativity (GR). 

In cosmology of homogeneous and isotropic Universe, the points in the space with 
coordinates r move as given by a dependence of scale factor on the time, a = a{t), so 
that velocity v and acceleration g of material point are determined by 

X = a{t)r, v = x = dr, g = x = dr, (1) 

where the dot over a symbol denotes the derivative with respect to time. Due to 
the homogeneity and isotropy, the evolution equations do not contain the co-moving 
coordinate r. Therefore, the evolution law is universal, and it can be considered at 
small r, i.e. in the region of applicability of the non-relativistic mechanics {v — > 0) and 
gravity law by Newton, namely, the acceleration of free falling along the radius vector 
g = dr reads off 

GMkt 



where we introduce the Komar-Tolman gravitational mas^il for a spherically symmetric 
distribution of the matter with the energy density p and pressure p inside the sphere 
with radius 

An 

A^XT = y (p + 3p)|a;|3, (3) 

so that 

d AnG 



. (p + 3p), (4) 

(X o 

that naturally does not include the co-moving coordinate r. 

In addition, the adiabatic law of energy conservation for the substance with the 
energy density p and pressure p reads off 

d£ = -pdV (5) 

so under the substitutions of £^ = pV(t), V{t) = a^Vo, we get 

p + 3H{p + p) = 0. (6) 

I The source of gravitational field in Einstein's equations is the tensor of 2{T^i, ~ g^u), where the 
symbol T corresponds to the tensor of energy-momentum and its trace, and g^^ denotes the metric. 
Therefore, the static gravitational potential is determined by not the energy density alone, but by the 
temporal component of the source. The integration of temporal component over the volume gives the 
Komar-Tolman mass, exactly. 
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Then, after multiplying by 2aa we can easily integrate and get the Friedmann 
equation 

/^^^^(P.f). (T) 

where we have introduced the Hubble constant H = a/ a, while in GR the constant 
of integration p^; is presented by the sum of contributions given by a constant spatial 
curvature and a matter with equation of state p = — | p- These terms do not contribute 
to the Komar-Tolman mass, i.e. they do not produce the gravitational force. 

Thus, we get the equations of GR for the evolution of homogeneous and isotropic 
Universe (|6]) and ([7]). Then, the description of such evolution is reduced to the non- 
relativistic dynamics at r — > 0. 

The straightforward application of MOND paradigm at r — )■ implies the transition 
to the regime of modified law of gravity, when 



9^9 

hence. 



9_ 

9o 



at 15-1 < 5-0, 



-W = ^ P + 3p — . 8) 

a 3 r 

Notice that a direct observation of the scale factor evolution by the registration of type 
la supernovae [3 El [9], [TOl [11] reliably shows that the accelerated Universe expansion 
(a > 0) at present has changed the decelerated expansion (a < 0), which took place 
at the redshiftl§] z > zt, where Zt ~ 0.4 — 1.0 is extracted in the procedure of fitting 
the deceleration parameter q{z) = —d/{aH^) by the linear function q{z) H- go + ^Qo- 
Consequently, in vicinity of transition point Zt, where a = 0, the validity of MOND in 
cosmology is quite justified. 

However, eq. (E]) at = const is inconsistent with the condition of homogeneity and 
isotropy of matter, since the actor ^ explicitly depends on the co-moving coordinate r, 
hence, such the modification of gravity law is inherently related to the inhomogeneity of 
matter distribution. In particular, eq. ([H]) results in the vacuum catastrophe: the vacuum 
with Pa = —PA = const is to be unstable, namely, due to eq. ([8]) the inhomogeneous 
distribution of matter is initiated. Thus, the straightforward transfer of MOND to 
the cosmology is theoretically inconsistent. Nevertheless, keeping in mind that MOND 
is related to the spatial inhomogeneity of matter distribution, we can easily suggest 
its cosmological extrapolation for the homogeneous and isotropic Universe, indeed: by 
setting 

^0 fl-o = fl-ola^l, (9) 

we arrive to the evolution equation at \d/a\ q'q 

d \d\ AttG , . , , . 

-— = -^(P + 3p)^? , (10) 
a a 3 

§ The redshift is related to the scale factor hy z ~ — 1. 
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which is consistent with the initial conditions imposed on the matter as well as with the 
vacuum stability. 

On the other hand, the visible large scale structure (LSS) is seen at the angle 
distances of ~ 1.5 — 2° ~ 'k)~ ^ radians, that corresponds to |a3|Lss ~ ^^lss ■ 
where the Hubble horizon ~ at the current value [B] 



km 

Iq ^ (69 - 72) 



Empirically in MOND [2] 
Ho 

271 



Mpc 



^o-f^, (11) 



so that if we put 



go ~ fi^ola^lLss, 



then we expect 



9o^7j-^Hl (12) 



It is convenient to parameterize the quantity Qq in the form 

9o = KoHl (13) 
whereas by the order of magnitude 

Ko ~ 10. (14) 

Thus, we can assume that the change of regime in the cosmological extrapolation 
of MOND as given by (j9]) is caused by the presence of spatial inhomogeneity of matter, 
that allows us theoretically to conform the existence of two regimes for the critical 
acceleration of gravity. 

In this paper we investigate the scheme proposed in and (1131) for the cosmological 
extrapolation of MOND with the parametrization in the transition region as derived in 
the holographic description of gravity as the entropic force [121 IE] at low temperatures 
(the MOND regime) [H]. This way essentially differs from the relativistic theory of 
gravity for the MOND paradigm offered by J.Bekenstein in [15], wherein he constructed 
the action satisfying a set of requirements, which include the appropriate limit to the 
MOND regime as well as, in principle, a correct description of cosmological effects, 
too. The modification is very specific due to additional gravitational scalar and vector 
fields in the action. Additional fields appear in the modified gravity formulated by 
J.W.Moffat in [16], too. Then, the gravitational constant becomes the field satisfying the 
equations of motion. So, J.W.Moffat has obtained the modification of gravity law [17] 
similar to the MOND, but it is distinguishable from the MOND, while some cosmological 
effects can be treated in agreement with observations [18j. Thus, approaches in [15] and 
[16] represent the axiomatic attempts to construct the consistent theory of modified 
gravity from the primary principles of field theory, and they are able to produce 
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the framework for complete calculations including the propagation of cosmological 
perturbations in the relativistic way, of course. We follow a more pragmatic way: 
without an introduction of any new additional notions, the modification of gravitational 
dynamics in the form of ([9]) and ( IT3|) can be consistent with the MOND at galactic scales, 
and it is enough for the investigation of modified cosmological evolution essentially 
different from the evolution in the framework of general relativity if we consider the 
ordinary components of matter. Therefore, the extrapolated MOND could exhibit some 
general features, which are common for meaningful extensions of general relativity as 
they tend to exclude the dark matter from the cosmology dynamically. In this way, 
we understand that the specific additional gravitational degrees of freedom defined 
explicitly as in |15j and [T6j should be inevitably considered if we try to take into 
account the modification of theory for the propagation of matter inhomogeneity in 
cosmology, hence, in the construction of full range description of large scale structure of 
Universe, that includes the matter power spectrum and anisotropy of cosmic microwave 
background radiation. In the general relativity, that is the dark matter, which is 
responsible for the actual description of all features in the large scale structure due 
to corrections to the matter distributions during the cosmological evolution, but it 
is commonly known that basic cosmological characteristics such as the accelerated- 
decelerated regimes of Universe expansion, and the visible angular scale of matter 
inhomogeneity in the sky are mainly given by the appropriate evolution rate of Universe, 
and we show that such the evolution can be achieved in the framework of cosmological 
extrapolation of MOND without any dark matter. This extrapolation does not provide 
us with the complete theory of inhomogeneity propagation (especially, because such 
the theory is expected to be nonlinear in contrast to the well-working ordinary linear 
perturbation theory in the general relativity), hence, we do not suppose to reproduce 
the full description of angular dependence of cosmic microwave background radiation 
subject to the given primary spectrum of matter distribution, of course, or the big- 
bang nucleosynthesis, for instance. Nevertheless, the baryonic matter inhomogeneities 
themselves can produce the potential gravitational wells for the further concentration 
of matter, though such the wells would be affected by the sound propagation in the 
photon-electron-baryon plasma, and this evolution of inhomogeneities will be mainly 
determined by the strong field regime of MOND, which is consistent with the general 
relativity, while the nonlinear effects are expected to be subleading. Thus, we can expect 
the correct reproduction of main features of large scale structure in the framework of 
cosmological extrapolation of MOND. So, our preliminary studies show that the first 
acoustic peak in the cosmic microwave background anisotropy is successfully obtained 
within the cosmological extrapolation of MOND under the very standard settings of 
primary spectrum of matter power without any dark matter, of course, while the 
second peak can be also correctly described by an appropriate variation of the matter 
power spectrum, but the subleading structures need the construction of full theory of 
cosmic perturbations beyond of purposes of our paper. So, the analysis of anisotropy 
in the cosmic microwave background radiation would be presented elsewhere, since it 
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requires a deep and comprehensive investigation in details. The similar note can be 
addressed to the problem of nucleosynthesis, too. Here, we present the general scheme 
for the cosmological extrapolation of MOND and its applications to the basic features 
of Universe evolution. 

In this paper, we apply the minimal scheme of ([9]) and (fT3l) to the description of data 
on type la-supernovae. We show that in the redshift region accessible for observations, 
such the astronomical data can be reliably described in the proposed model if there 
are the cosmological constant and baryonic matter, only, without any dark matter, 
whereas the numerical estimate of ( IT^ is correct. While comparing with the Friedmann 
cosmology at presence of dark matter, we can evaluate the imitated contribution of dark 
matter pm to the critical density of pc = 3Hq/8ttG. The ratio of = Pm/Pc to the 
baryonic contribution Qb = Pb/Pc is roughly determined by Kq in order of magnitude. In 
accordance to the standard Friedmann cosmology-model with the cold dark matter and 
cosmological constant (ACDM), the parameters correspond to ^ 0.27 at fib = 0.045. 

Further, we analyze the large scale structure of inhomogeneity in the matter 
distribution of Universe in calculations of baryon acoustic oscillations, measured by 
observing the luminous matter in the sky versus the redshift [I9], and in estimation of 
"acoustic scale" in the CMBR anisotropy [6]. The estimates obtained in the framework 
of cosmological extrapolation of MOND agree with the measured values of LSS without 
any postulation of hypothetic dark matter: it is enough to adjust the value of baryonic 
matter in the presence of cosmological constant. Then, the density of baryonic matter 
is roughly twice its value in ACDM of GR. 

After modeling the acoustic scale in the CMBR anisotropy we calculate the critical 
acceleration qq in the framework of the cosmological extrapolation of MOND and find 
preferable values of cosmological parameters of the Universe (the density of baryonic 
matter). We analyze the evolution into the future, too. 

In Conclusion we summarize our results and discuss its meaning. 

2. MOND parameters from data on supernovae 

At present, measuring the stellar magnitude of type la supernovae versus the redshift 
gives the only direct observation of scale factor versus time in cosmology. Such the 
observation does not base on any model of Universe expansion, so that even the theory 
beyond the Friedmann equation can be confronted to the experiment without any 
additional assumptions. In this respect, such the data are unique and they are suitable 
for testing the cosmological extrapolation of MOND. 

The stellar magnitude p depends on the redshift according to the formula 

p = Pabs + 5 log^Qdiiz) + 25, (15) 

where pabs is the absolute stellar magnitude, i.e. the stellar magnitude of light source 
at the distance of 10 pc. The photometric distance di measured in Mpc, is determined 
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by the Hubble constant evolution 

z 

dL{z) = {l + z) 1 (16) 



where c is the speed of light. In [7] the following dependence of the deceleration 
parameter on the redshift has been used as the working hypothesis: 

q{z) ^ qiin. = qo + zq'o, (17) 

that allows it to describe the supernovae data at values 

go = -0.85 ±0.35, qQ = l.8±l, (18) 

essentially different from zeros, whereas the parameters significantly correlate (see details 
in the original paper [7j). 

According to (fT7|) . the Hubble constant 

Hnn. = Ho{l + z)i+^»-^o exp{q',z}, (19) 

that should be reasonably compared with the case of zero deceleration parameter, when 

H^,i = Ho{l + z). (20) 

From eqs. (fT5|) - (l20|) we see that normalizing the data at low redshifts, i.e. measuring 
the value of Hq, allows us to plot a fiducial dependence of supernovae stellar magnitude 
Ho with the deceleration parameter equal to zero, and further, systematically to study 
the deviations from this dependence versus the increasing redshift. 

Data from [7] are shown in figures [2] and [3l They confidently signalize on the 
accelerated expansion of Universe at present as well as on the existence of transition 
from the decelerated expansion to the accelerated one aX z = zt ^ 0.5. The best fit of 
data takes place at parameter values listen in f lTSj) . as it shown in figure [3] by the dashed 
curve. 

In the cosmological extrapolation of MOND we use the following equation for the 
acceleration at the presence of both the cosmological term with the energy density p\ 
and baryonic matter with the density pb- 

-^(iTTt] =^(2PA-P.), (21) 
a \ o |a| / 6 

where the region of transition from the Friedmann regime to the MOND is described 
by the 1-dimensional Debye function as accepted in [H], hence, 

DW = i/^^. (22) 

X J exp{y} - 1 



so that at the acceleration greater than critical one, i.e. in the limit of standard 
cosmology at |d|/a ^ (^q, we get 

D a 
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Figure 2. The stellar magnitude /i of type la supernovae versus the redshift z: the 
bottom dotted curve presents the case of zero deceleration parameter, the solid curve 
shows the result of cosmological model of MOND, the top dotted curve depicts the fit 
in ACDM. The supernovae discovered by the ground based telescopes are marked by 
rhombuses, the circles show the supernovae discovered by the Hubble Space Telescope 
(data from [7]). 



while in the deep regime of MOND, i.e. at \a\/a ^ q'q, asymptotically we find 

6 \a\ J ag'o 



and we arrive to eq. ( ITOl) . 

The critical density defined by 



n - ^ 



allows us to write down the baryonic density as 

fife 

Pb = Pc^, 

while the requirement of vanishing the acceleration at redshift Zt leads to the following 
expression for the cosmological contribution: 

Pa = ^pA(l + zt)^. 
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Figure 3. The difference of stellar magnitudes A(/x — /iq) for the type la supernovae 
versus the redshift z after subtracting the fiducial magnitude at zero deceleration 
parameter. The solid line refers to the cosmological MOND, the dashed line shows 
the best fit at the deceleration parameter q{z) — qq + zqq [7]. The bottom panel 
serves for the systematical illustration of supernovae data as averaged for the similar 
redshifts. The dotted and dashed-dotted lines represent the evolution with appropriate 
constant positive and negative deceleration parameter, respectively. 



Thus, the model for the cosmological extrapolation of MOND is described by equation 
combined with initial data 

Figures [2] and [3] illustrate the results of cosmological MOND at the typical value of 
Qb = 0.045 as in ACDM, while 

go = -0.775, zt = 0.375, ^ Kq = 16.7, (24) 

so that in the region of available data, the calculated profile is practically 
indistinguishable from the fit obtained at the linear dependence of deceleration 
parameter on the redshift as performed in [7]. Note, that the difference of stellar 
magnitudes presented in figure [3] is independent of the normalization of Hubble constant 
Hq, while data shown in figure [2] allows us to fix both the value of Hq at small z and 
the fiducial curve at a = 0. 

Notice that at present, the regime of deep MOND takes place in cosmology, so that 
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from eq. ( 12^ we find 

ql = ]^K^^l,{{l + Ztf-l). (25) 

Hence, we can see that the value of baryonic matter density, i.e the parameter of can 
get quite arbitrary variations in this test on the supernovae data, because the region of 
MOND apphcability near a = is completely described by the factor of K^VLb, entering 
in eq. f l25|) . while the Hubble constant at low redshifts is mainly determined by initial 
data, namely, by Hq and go, so that the influence of Vt^ value on the quality of data 
description is quite weak. For instance, the twice increase of baryonic fraction is almost 
invisible on the curve derived from the cosmological MOND as shown in figure [3l i.e. 
the region of low redshifts is weakly sensitive to fife. 

In the framework of ACDM the deceleration parameter and redshift of transition 
from the deceleration to the acceleration are given by the fractions of densities for the 
vacuum Vt^ and matter VLmi 

go = ^ - fiA, (1 + ztf = 2^, (26) 

so that in the case of spatially fiat Universe, when + = 1, "we get 

?o = ^fiM-l, {l + ztf = 2^—^. (27) 

The quality of fitting the supernovae data for the fiat Universe in ACDM at 
~ 0.73, VLm ~ 0.27 is slightly worse than in the case of the linear approximation 
( IT71) for the deceleration parameter, but it has got a little bit of better quality at 
^A ^ ~ 0.95, ^ (Im ^ 0.45 [7]. 

The parameters brought from ACDM into MOND give 

— = -— = ^28) 
fife go fiA — I fiM 

At I go I ~ 1 we find fi^f/fife ~ Kq. However, by taking into account the fact that 
the point of maximal confidence level of data fit in ACDM is shifted from the standard 
values fiA ~ 0.73, Qm ~ 0.27 to fiA •-)■ fiA ~ 0.95, Qm ^ ~ 0.45, hence, numerically 
fi^w/fife ~ 10, we deduce Kq ^ 7.5, that is in agreement with our estimates obtained in 
consideration of changing the regime of homogeneous matter distribution to the limit 
of inhomogeneity in MOND. 

Thus, in the framework of cosmological extrapolation of MOND we can reliably 
describe the type la supernovae data by making use of baryonic matter and cosmological 
constant, only, without any postulating the dark matter. In this way, the parameters of 
extrapolation is naturally conformed with the Milgrom's critical acceleration in MOND, 
applicable at scales, when the spatial inhomogeneity of baryonic matter distribution 
becomes significant. The cosmological extrapolation of MOND allows us, in fact, 
to calculate the present ratio of density of hypothetical dark matter to the baryonic 
density fiA^/fife in terms of the slope of critical acceleration in its dependence on the 
distance, i.e. the parameter of g^, expressed in units of Hubble constant at present. 
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Table 1. The comparison of observed parameters for the large scale structure of 
Universe with estimates in the framework of the cosmological extrapolation of MOND 
(the parameters are explained in the text). 



quantity 


rs{zd) 


Dv(0.35) 


Ia 


exp.[6l[T9] 


0.1905 ±0.0061 


0.1097 ±0.0036 


302.69 ±0.76 


MOND 


0.1924 


0.1149 


302.5 



Therefore, the supernovae data lead to the necessary introduction of dark matter if 
we keep the Friedmann evolution, i.e. if the evolution is determined by the Newton 
law of gravity, indeed, while the modification of gravity law at low accelerations in 
the case of inhomogeneous distribution of matter in accordance with MOND as well as 
for the homogeneous distribution of matter in cosmological MOND, gives the confident 
description of data without any dark matter. 



3. The large scale structure of Universe 

The Universe evolution essentially transforms the spatial distribution of matter, that is 
observed as baryonic acoustic oscillations (BAO) [12] and CMBR anisotropy ^. 

So, the results of survey of the matter distribution on the celestial sphere versus 
the redshift [19] are particularly reduced to the ratio of sound horizon for the waves 
in the baryon-electron-photon substance as written in co-moving coordinates 

t{z) 

r,{z) = ! Csdt (29) 







to the effective distancqj]] defined by 

Dv{z) = i^{l + zfD\{z)^y, (30) 
where is the sound speed in the medium, while the angular distance is given by 



The measurements give the value of rs{zd) / Dy{z) aX z = 0.2 and z = 0.35, which are 
shown in table [H Here Za is the redshift of epoch, when the interaction of baryonic and 
photonic components of substance via the Compton scattering of electrons off photons 
and Coulomb attraction of electrons to protons becomes negligible, i.e. when baryons 
stop to drag photons (the finish of drag epoch). 

II We consider the flat space. 
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The speed of sound is given by the following expression: 

c 1 



(32) 



73 yrr^' 

where R is the baryon-photon ratio depending on the redshift 

4 4 (1 + z)VL^ ^ ' 

Introducing the redshift z^q for the moment, when the density of non-relativistic matter 
Pm equals to the density of radiation (photons and neutrinos), allows us to write down 
the formula for the sound horizon as 

r,(.) = A /Zin^i±5±V^?±5^, (34) 

keq y Req 1 + a/ Req 

where 

Analytical parameterizations of numerical calculations for Zf.q, Zd and R are given in |20] . 
wherein the dependencies on the CMBR temperature T = ■ 2.7 K and parameters 
VLi)h^, VlMh"^ with h being the coefficient for the Hubble constant written down as 
Hq = h- 100 km ■ s~^/Mpc, are explicitly presented. Thus, the Universe evolution at low 
redshifts determines the angular distance -D^ in eq. (13T]) . which depends on the initial 
data on Hq and qo, while the sound horizon and redshift Zd are determined by fractions 
of baryonic and dark matters in the energy budget of Universe. Since the evolution at low 
redshifts in the cosmological extrapolation of MOND conforms with the supernovae data 
in the same way as it was done in ACDM, the angular distances are almost coincident 
in both these models, hence, the description of BAO data is determined by fitting the 
sound horizon r^. In the case of the cosmological extrapolation of MOND all of non- 
relativistic matter is identified to baryons, i.e. Qm = ^b- Hence, we find Qhh"^ ^ 0.057, 
i.e. it is twice and a half greater than the baryonic fraction in ACDM. At such the value 
of baryonic density we get the agreement with the BAO data (see table [I]). In numerical 
estimations we have put h = 0.71, T = 2.725 K and go = —0.775, Zf = 0.375 as it was 
in the previous section. Then, Qi, = 0.113. Such the doubling can point to necessity of 
accounting for the baryonic matter outside luminous stars and visible gas, that is also 
numerically consistent with the lack of visible baryonic matter in some galactic clusters, 
which dynamics is described in MOND without any dark matter (see [21 13]). 
At the same parameters we estimate the "acoustic scale" 

Ia = {1 + z^} -— , (35) 

which is measured in the spectra of CMBR temperature anisotropy by locations of peaks 
of temperature fiuctuations versus the multipole number in WMAP data [6]. Here Zj, 
is the redshift of decoupling, when due to the recombination of electrons with protons 
the medium becomes transparent for photons (see analytical approximations for z^ in 
terms of baryonic density, matter density and Hubble constant in t21j-)- Calculating the 
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value of Da{z^)^ we take into account for the contribution of relativistic particles into 
the Hubble constant, since 

l + z 

Pt + Pu = Pbr— . 

As we can see from table [H the cosmological extrapolation of MOND is safely able to 
describe the large scale structure data given by the measurements of both the baryon 
acoustic oscillations of visible matter and the anisotropy of CMBR0. 

It is interesting that our argumentation about changing the MOND regime at 
distances of the order of the scale for the spatial inhomogeneity of matter distribution, 
in fact, is confirmed by the calculation of acoustic scale which points to the angular 
size of inhomogeneity about 360°//^ ~ 1.2° in consistency with our initial assumptions. 

4. The Milgrom's acceleration 

Since the critical Milgrom's acceleration in MOND is experimentally measured with the 
high accuracy in the gas-reach galaxies |5] 

go = (1.21 ± 0.14) ■ 10"^°m/s^ (36) 

it makes sense to use this information in order to refine of parameters in the cosmological 
extrapolation of MOND. 

Indeed, within our approach, the Milgrom's acceleration is determined by the 
angular scale of inhomogeneity in the matter distribution 66, so that in Hubble units 

4^ = Ko69. (37) 

However, the acoustic scale gives 

Se = ^. (38) 

Then, by expressing Kq in term of deceleration parameter go and redshift Zt for the 
transition from the Universe deceleration to its acceleration in accordance with fl2^ . we 
find 

9o _ 477 ql 



(39) 



% We have not considered the "shift parameter" defined in [22j, !R = y^ilMH^{l+z^)DA{Zf:)/c, since, as 
authors of 22 have pointed out, this quantity can not be directly extracted from the angular spectrum 
of CMBR anisotropy because of a poor accuracy of measurements at low values of the multipole 
number, i.e. at large angles of correlations, wherein the shift parameter makes a significant infiuence 
on the amplitude of spectrum, while the high multipoles are not sensitive to Ji. It means that the 
choice of model for fitting the spectrum in the region of acoustic oscillations a.t I ^ Ia determines the 
extrapolation to low values of multipoles, but it is not critical for the data at low multipoles because of 
poor accuracy of data at low I <^Ia- Consequently, the value of shift parameter given in [6l is actually 
the model dependent-extrapolation, indeed. 
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that allows us to plot the dependence of Milgrom's acceleration on the baryon density 
at fixed values of go aiid Zt, which are determined by the supernovae data. Further, we 
can compare the result with the empirical value of 

4^ = 0.180 ± 0.022, (40) 

C-Ho 

known with the 12%-accuracy after taking into account for the uncertainty of the Hubble 
constant. 
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Figure 4. The Milgrom's acceleration go/cHo calculated in accordance to 
at the deceleration parameter go = —0.853 and the redshift of zero acceleration 
Zt — {0.338, 0.375, 0.413} (the curves from top to bottom, respectively) as well as 
at given acoustic scale Ia = 302.6 (the dashed line) versus the baryon fraction ili, in 
comparison with the experimental value (the shaded band). 



The comparison is shown in figure HI It turns out that the agreement of data 
on the supernovae (SN), baryonic acoustic oscillations (BAO), acoustic scale (WMAP) 
and Hubble constant (Ho) with the Milgrom's acceleration is reached at the deceleration 
parameter go ~ —0.853 and accessible 10%- variation of the transition redshift zt ~ 0.375. 
In this way, the account for the small error in the value of acoustic scale leads to 

Qb = 0.115 ±0.012, 

go = (1.10 ± 0.03) ■ lO^^V/s^, 

(41) 



DviO.2) 



0.185 ±0.006, 
0.112 ±0.004. 



Dv-(0.35) 

Thus, we successfully conform the logical consistency of two MOND regimes: the 
cosmological limit and the locally inhomogeneous case; and we extract the confident 
intervals for the Universe parameters. 



Cosmological extrapolation of MOND 
5. The evolution into the future 



16 



In the future a{t) ^ 1, and the equation of cosmological MOND is reduced to 

1 



a \a\ 
a a 2 
hence, according to ( 125) 

a \a 
a a 



(42) 



(43) 



l + Ztf-V 

so that the limit is the de Sitter Unverse with a{t) = a-^exp[H\(t — t^)], where the 
cosmological Hubble constant 



Since Zt > and Ha ^ Hq, 



% < 1, 



1. 



(44) 



(45) 



Therefore, the future epochs at Ha ^ Hq are constrained by condition ( HSil as illustrated 
in figure O 



Zt 

0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
01 



0.0 
0.60 



Ha < Hq 








l --^'^^^^^^^^^ Ha > Hq 





m 



Figure 5. Delimiting the future epochs versus the deceleration parameter qq and the 
redshift of zero acceleration zt. The preferable fit of WMAP+BAO+SN+_ffo+5o data 
is depicted by the crossing of the shaded bands. 

Substituting of go = —0.853 as obtained by fitting the data, into eq. ( H5|l gives 

Zt ^ 0.54, 

so that according to our estimates for Zt in the previous section, it means that to the 
end of evolution Ha > Hq (see figure [5]). 

It is interesting to compare the cosmological term evaluated in MOND in terms of 
Ha with its value in GR: 

Hl=^-^PA = ln,H^oil + Ztr, (46) 
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so that the equahty of //a = H\ takes place at 

o2 4gg 1 
" (l + z,)' (l + ^,)3_i- 

Numerically 

fib < fife. 

Therefore, the limit of de Sitter Universe and, hence, the problem of small-scale 
cosmological constant cannot get an adequate description in GR, while in order to 
consider the vacuum we need essentially to modify the equation for the connection 
between the gravity and matter. 

6. Discussion and conclusion 

In this paper we have conformed the MOND to its cosmological extrapolation, 
that allows us to eliminate hypothetical dark matter not only at scales of spatial 
inhomogeneity of baryon matter distribution in galaxies and galactic clusters, but also 
in the description of Universe evolution up to the redshift of 2; ~ 1.8 in the observation 
of type la supernovae. In this conformation, the angular size of inhomogeneity in the 
large scale structure of Universe is essential. 

Another question, related to the cosmology and modification of gravity, is the 
necessity to postulate the dark matter in order to describe the CMBR anisotropy in 
ACDM. To our point of view, this problem is also closely connected to the usage 
of Friedmann model of expansion at redshifts from zero to 2;*, when the CMBR 
inhomogeneity was formed. Clearly, parameters such as the contribution of cosmological 
term into the energy budget of Universe and the fraction of dark matter can be mostly 
significant at low redshifts, exactly when the modification of gravitational law becomes 
remarkable at low accelerations of expansion. However, the forming of CMBR essentially 
depends on the proportion between the baryons and dark matter. On the other hand, 
the evolution of the baryonic density versus the redshift is completely determined by 
the conservation law of matter, hence, the conformation of the baryon densities in 
both ACDM and MOND allows us to generate similar conditions for the interaction 
between photons and electrons, which density correlates with baryons. The problem 
of considering the anisotropy evolution versus the redshift at the whole interval in 
the framework of cosmological MOND expects a solution, i.e. the question about a 
detailed investigation of CMBR anisotropy still remains open within the cosmological 
extrapolation of MOND. In this respect we have to mention the field theory-approaches 
in [15] and [I6], which can, in principle, provide us with complete calculations for 
the propagation of perturbations, whereas the theory by J.Bekenstein [15] includes the 
MOND regime for inhomogeneous distribution of matter. However, we follow another 
motivation based on the minimal extension of MOND that satisfies the reasonable 
constraint on the transition to the cosmology of homogeneous distribution of matter. 
In this way, we note that some integral quantities, which values are deduced by the 
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consideration of CMBR anisotropy, are derived under the assumption of Friedmann 
evolution and they cannot be straightforwardly transferred to the modified law of gravity. 
This note is also valid as concerns for the observation of large scale structure, i.e. in 
studying the baryon acoustic oscillation^ (BAO) caused by the spatial inhomogeneity 
^19j . Nevertheless, we have analyzed the "angular scale" in the CMBR anisotropy and 
parameters of baryon acoustic oscillations in the framework of cosmological MOND and 
showed that MOND can confidently describe these cosmological phenomena. Then, we 
can conclude that the cosmological MOND has exhibited its positive potential for the 
adequate description of Universe evolution at low redshifts as well as up to the redshifts 
of CMBR forming. We have to note, of course, that the purpose of this paper has 
been not a complex fitting the type la supernovae data, baryon acoustic oscillation and 
CMBR anisotropy, but the aim has been the demonstration of opportunity to get the 
successful application of cosmological MOND to those problems. 

Finally, we have compared the critical acceleration calculated in the framework of 
cosmological MOND with its empirical value and extract the intervals of cosmological 
parameters, wherein there is the good agreement of the model with observations of 
WMAP+BAO+SN+ Hq+Qo- In addition, we have found that the modification of gravity 
law is essential in studying the problem of cosmological constant. 

Thus, the preliminary estimates performed in the framework of cosmological MOND 
conceptually permit for the confidential description of Universe evolution and point 
to the necessity to work out a cumbersome analysis of complete base of the data on 
the supernovae, large scale structure, baryogenesis and abundance of elements in the 
Universe etc., for the accurate extraction of cosmological parameters, in particularly, 
the baryonic fraction of energy. Nevertheless, in this paper we have shown that 
hypothetical dark matter can be completely excluded from the cosmology by the 
appropriate modification of gravity law at accelerations less than the critical value in 
analogy to the case, when it was successfully done at the galactic scale in MOND. 

+ Notice, in the framework of ACDM the spatial inhomogeneity of both the matter and dark matter 
are generated, say, by the primordial spectrum of inhomogeneity as given by quantum fluctuations of a 
scalar filed having caused the Universe inflation (23l [24l [25l [26l [27] . Therefore, the baryonic and dark 
matter initially has got common locations of compression-decompression regions, that is consistent 
with the assumption about the joint concentration of baryonic and dark matter. However, the baryon 
oscillations having the electromagnetic nature, suggest that the inhomogeneities of the photon-electron- 
baryon medium propagate as sound waves decoupled from the dark matter (except the interaction due 
to gravitational forces), while the inhomogeneities of dark matter do not propagate at all, since there 
are no sound in the dark matter. During the evolution after the baryon acoustic oscillations stopped, 
the gravitational interaction attracts the centers of baryonic and dark matter concentration, that results 
in a partial nearing the locations of matter lumps. The transfer function of initial perturbations in 
the energy densities consists of two different parts given by the terms of baryonic and dark matter. 
Consequently, initially coinciding regions of concentration for the baryons and dark matter become 
spatially separated, that creates different centers of gravitational contraction of baryons and dark 
matter in the future evolution. Thus, the existence of baryon acoustic oscillations challenges the dark 
matter hypothesis by asking for a natural explanation for the coincident locations of baryons and dark 
matter in galaxies, that rises the additional problem in a dynamical derivation of the Tully-Fisher law 
in the framework of standard cosmology. 



Cosmological extrapolation of MOND 



19 



Notice, there are several papers, wherein the influence of modified law of gravity 
to the Universe evolution at low accelerations was considered [281 [29], however, in those 
papers the acceleration is evaluated by the temperature of apparent horizon Th = 11/211 
giving gii = 2'kTh = H, while the critical acceleration is the constant value independent 
of distance. In [29] this approach follows the MOND at the linear dependence of factor 
modifying the acceleration dr on ratio g/gn ^ 1- Hence, it is clear that such the 
modification of gravity law at the Hubble horizon is extrapolated inside the Hubble 
sphere, that can be never conformed with our consideration. 
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